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ABSTRACT 


The  magnitude  and  characteristic  of  the  effect  of  cloud-droplet 
microwave  dispersion  on  line  integral  ref ractometer  (LIR)  measurements 
was  determined  for  environmental  conditions  measured  at  the  Eastern  Test 
Range  on  August  7,  8,  and  9,  1963. 

It  was  definitely  determined  that  the  presence  of  cloud  droplets 
would  introduce  errors  in  the  LIR  measurements,  but  that  on  each  of  the 
above  days,  there  were  regions  of  the  sky  where  the  error  introduced 
was  small  compared  to  the  total  refraction  correction. 

The  report  describes  technical  areas  where  basic  data  were  inade¬ 
quate  and  discusses  the  impact  of  assumptions  made  in  these  areas  on 
the  calculated  values. 

Basic  work  in  this  area  at  The  MITRE  Corporation  in  1962  is  de¬ 
scribed  in  Appendix  I. 
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SECTION  I 


INTRODUCTION 

The  National  Academy  of  Science  Ad  Hoc  Panel  on  Electromagnetic 
Propagation  reported  in  1962  that  atmospheric- ref raction  effects  were 
one  of  the  major  limitations  of  missi le- range  metric-measurement  capa- 
bil ity . 

In  1963  work  was  initiated  at  The  MITRE  Corporation  on  a  new  ap¬ 
proach  to  solve  this  problem,  the  line  integral  ref ractometer  (LIR) 
technique.  This  approach  utilizes  the  differential  phase  shift  between 
two  coherent  signals,  15.6  and  31.2  gHz,  as  a  measure  of  the  integral 
of  the  water  vapor  component  of  the  refractive  index  over  the  actual 
transmission  path.  Two  other  coherent  signals,  45  and  90  gHz,  are 
used  in  a  similar  manner  to  measure  the  effects  of  the  dry  constituents 
of  the  atmosphere. 

Basic  theoretical  studies  indicated  that  water  droplets  in  clouds 
would  induce  error  in  the  LIR  measurements;  however,  field  tests  were 
undertaken  to  determine  the  feasibility  of  utilizing  this  technique  to 
meet  current  and  future  Air  Force  needs. 

Some  experimental  data  concerning  the  magnitude  of  the  water-droplet 
effect  on  the  differential  phase  measurements  at  15.6  and  31.2  gHz  were 
obtained  during  a  field  test  of  the  LIR  at  Lake  Winnepesaukee ,  New 
Hampshire.  The  amount  of  applicable  data  was  limited  because  the 
occurrence,  frequency,  and  the  type  of  clouds  in  the  fixed  measurement 
path  was  one  of  the  uncontrollable  variables  encountered  during  the 
relatively  short  testing  period.  While  the  test  results  confirmed  that 
the  effect  predicted  by  theory  was  indeed  present,  the  results  did  not 
provide  much  information  concerning  the  magnitude  of  the  effect  to  be 
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expected  for  an  LIR  measurement  system  at  the  Air  Force  Eastern  Test 
Range  ( ETR)  . 

No  attempt  was  made  to  determine  how  the  cloud  structure  prevailing 
on  the  3  days  in  August,  1963,  compares  with  that  prevailing  during 
other  times  of  the  year;  while  such  data  may  be  available,  it  was  con¬ 
sidered  to  be  beyond  the  scope  of  this  study.  However,  Dr.  Robert 
Cunningham  (AFCRL)  has  indicated  that  a  cloud  pattern  of  scattered 
small  cumulus  clouds  over  the  land  area  and  only  a  few  small  cumulus 
clouds  (or  none  at  all)  over  the  water  is  typical  for  a  very  high  per¬ 
centage  of  the  days  in  the  Florida  ETR  area.  A  photograph  (Figure  l), 
taken  from  Gemini  in  July,  1963,  shows  a  similar  cloud  pattern. 

A  brief  summary  of  the  LIR  feasibility  test-result  conclusions  is 
included  from  the  final  report  L .  The  test  program  demonstrated  that 
it  is  feasible  to  use  the  microwave  dispersivity  at  13.6  and  31.2  gHz 
as  a  precise  measure  of  the  integrated  refractivity  of  the  more  variable 
component  of  the  atmosphere,  water  vapor.  Measurements  were  also  made 
at  45  and  90  gHz  to  obtain  information  concerning  the  oxygen  in  the 
atmosphere. c ^ 

The  2  ppm  accuracy  of  the  LIR  technique  represents  a  significant 
improvement  over  other  measurement  methods,  and  data  analysis  indicates 
that  accuracy  objectives  can  be  achieved  with  a  smoothing  time  of  10  to 
1000  seconds  (with  a  most  probable  time  of  100  seconds).  The  100-second 
period  of  time  corresponds  to  a  spatial  resolution  of  200  meters  at  a 
typical  crosspath  veloctiy  of  2  meters  per  second.  In  general,  theory 
indicates  that  the  integrating  time  will  decrease  as  the  angular  rate 
of  a  sweeping  beam  increases. 

The  field  test  measurements  demonstrated  that  the  LIR  technique 
could  be  utilized  to  achieve: 
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Figure  1.  Cloud  Distribution  at  ETR  (1965) 


(a)  a  2-  to  10-fold  improvement  in  accuracy  over  current 
state  of  the  art, 

(b)  a  20-  to  1000-fold  response- t ime  improvement  over  current 
state  of  the  art,  and 

(c)  valid  measurements  in  a  turbulent  atmosphere. 

Theory  exists  to  extend  the  fixed  path  results  to  a  moving  path  but 
experimental  verification  is  necessary.  The  LIR  instrumentation  used 
in  the  tests  was  designed  for  these  purposes  only  and  was  more  elaborate 
than  needed  for  an  operational  system;  however,  equipment  modification 
is  possible  to  yield  a  further  improvement  in  accuracy  and  a  reduction 
in  size  and  weight.  While  the  results  of  the  fixed  path  field  test 
measurements  prove  the  feasibility  of  the  basic  technique,  they  do  not 
afford  a  confirmation  of  the  rotating/sweeping  beam  effects. 
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SECTION  II 


BASIC  ASSUMPTIONS  AND  DATA 


The  magnitude  of  the  droplet  effect  is  dependent  on  the  density, 
temperature,  size,  and  distribution  of  the  droplets. 

The  calculations  are  based  on  the  assumption  that  all  the  liquid 
water  measured  was  due  to  cloud  droplets  and  not  rain  or  ice.  If  the 
liquid  water  had  been  in  the  form  of  rain,  the  effect  would  have  been 
much  more  detrimental;  if  it  had  been  in  the  form  of  ice,  the  effect 
would  have  been  less  detrimental.  In  addition,  the  dimensions  of  the 
droplets  was  assumed  to  be  small  compared  to  the  wavelength  of  the 
signal.  The  data  used  was  the  only  information  available  with  the 
detail  necessary  for  this  study.  The  visible  cloud  cross  sections 
for  the  3  days  are  shown  in  Figures  2,  3,  and  4. 


The  value  of  the  refractive  index,  n,  of  a  cloud  of  droplets  in 
the  atmosphere  as  a  function  of  percent  volume  content  of  water,  may 
be  expressed  as 


n 


*  *  *4 


-  1 


nd  +  2 


x  10 


(1) 


where 

M  =  the  equivalent  water  density  of  the  atmospheric  medium  in 

,  3 

grams/cm  , 

3 

P  =  the  density  of  liquid  water  in  grams/cm  ,  and 

n  =  the  refractive  index  of  the  water  molecule, 
a 


The  sensitivity  of  the  refractive-index  value  for  normal 

3 

about  1.5  refractivity  (N)  units  per  1  microgram  per  cm 
of  equivalent  water  density.  The  value  of  n^  for  normal 
on  the  order  of  7. 


conditions  is 
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Figure  2.  Cloud  Cross  Sections  for  August  7,  1963,  Valkaria,  Florida 
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Figure  3.  Cloud  Cross  Sections  for  August  8,  1963,  Valkaria,  Florida 
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Figure  4.  Cloud  Cross  Section  for  August  9,  1963,  Valkaria,  Florida 


The  complex  index  of  refraction  consists  of  dispersive  (real)  and 
attenuative  (imaginary)  components,  as  shown  in  Equation  (2). 

Refractive  Index  =  dispersive  +  j  attenuative  (2) 

or 

n  =  m  +  jk 

It  is  evident  from  the  sample  values  listed  in  Table  I  that  the 
characteristics  are  a  function  of  both  temperature  and  wavelength. 


The  variation  of  the  real  part  of  the  dielectric  function,  as  a 
function  of  temperature  and  frequency,  is  illustrated  graphically  in 
Appendix  I,  Figure  18.  The  two  frequencies  of  particular  interest  for 
the  LIR  application  are  15.6  gHz  (1.8  cm)  and  31.2  gHz  (0.9  cm).  The 
dispersive  effects  of  cloud  droplets  on  LIR  measurements  can  be  estimated 
from  the  difference  between  the  real  parts  of  the  refractive  index  at 
the  two  frequencies. 

Dispersion  =  1,5  x  10 


pf  116 


(fld2 


i) 


(n  2  +  2) 
d 


15.6  gHz 


(3) 


-  Re 


(n,  -  1) 

d 


( n . 


2) 


J31.2  gHz 


where 


m 


i  k. 


The  term 


Re 


( n  2  -  1) 


can  be  expressed  in  terms  of  the  real  and  imaginary  value  of  the  re¬ 
fractive  index  of  water, 
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TABLE  I. 


Values  of  the 
Literature  ^ 


Characteristics  of  Water  and  Ice  Crystals  Found  in  the 

3,4 


Temp.  (°C) 

\  (cm) 

m 

k 

Water 

-8 

.62 

3. 10 

1.77 

-8 

1.24 

4.15 

2.55 

-8 

3.21 

8.14 

2.00 

0 

.62 

3.45 

2.04 

3.49 

1.92 

0 

1.24 

4.75 

2.77 

0 

3.21 

7.80 

2.44 

10 

.62 

3.94 

2.37 

10 

.62 

4.08 

2.34 

10 

1.24 

5.45 

2.90 

20 

.62 

4.44 

2.59 

20 

1.24 

6.15 

2.86 

18 

.  9 

5.55 

2.85 

18 

1.25 

6.41 

2.86 

Ice  crystals 

0 

1.78 

0.002 

10 


Re 


1) 


2  2  2 
k  +  2)  +  4m  K 


(4) 


+  2) 


2  2  2 
(m  -  k  -  1)  (m  - 

2  2  2  2  2 
(m  -  k  +  2)  +  4m  K 


The  variation  of  dispersion  with  temperature  curve  in  Figure  5  is 
based  on  the  curve  in  Appendix  I,  Figure  20;  it  is  shown  as  a  dashed 
line,  whereas  the  straight  line  approximation  of  the  function  used  in 
the  computation  for  this  report  is  shown  as  a  solid  line.  There  was 
insufficient  data  available  concerning  the  temperature  regions  outside 
the  0-  to  10-degree  G  range.  It  was  assumed  that  the  curve  in  the  low- 
temperature  region  had  the  same  equation  as  the  linear  approximation  of 
the  function  in  the  0-  to  10-degree  C  range  and  that  it  has  a  constant 
value  of  0.01  in  the  temperature  region  above  10  degrees  C. 

N(t)  =  0.0310  -  (0.0021  x  A’Tn),  t  <  0°C,  (5) 

where 

ATn  =  °C  from  0°C, 

N(t)  =  0.01  t  >  10°C.  (6) 


The  effective  cloud  signal  at  elevation  angle  a  is  defined  as  D 
K  & 

D  =  N(t)  x  (MLC)  x  A£  x  360  x  10'B, 

a  S  n  n  n 


(7) 


n  =  1 


where 

N(t)  =  dispersion  at  the  measured  temperature  according  to 

the  graph  in  Figure  5  and  the  temperature  as  determined 
from  the  Rawinsonde  data  given  in  Appendix  II, 

3 

(MLC)  =  average  measured  liquid  content  (gm/m  )  for  the  section 
n  of  path  between  the  500-foot  grid  lines, 

£jL  =  distance  between  grid  line  n  and  grid  line  (n  +  1), 

K  =  uppermost  layer  containing  liquid  water  as  determined  by 

the  ETR  measurements. 
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Figure  5.  Dispersion  vs.  Temperature 


Measurement  data  on  the  amount  and  distribution  of  liquid  water  in 
a  cross  section  of  the  atmosphere  along  the  20-mile  east-west  baseline 
of  the  MISTRAM  system  at  Valkaria,  Florida,  and  extending  to  about  28,000 
feet  altitude  was  available  for  the  dates  of  August  7,  8,  and  9,  1963, 
from  measurements  made  by  Dr.  R.  L.  Cunningham  (AFCRL);  a  detailed 
description  of  the  measurement  program  is  given  in  References  5  and  6. 

The  temperature  vs.  altitude  profile  was  obtained  from  Rawinsonde 
measurements  up  to  20,000  feet  made  at  the  MISTRAM  Central  (MC)  site. 
These  are  tabulated  in  Appendix  II  together  with  the  liquid-water 
measurement  data.  Temperature  data  above  20,000  feet  was  obtained  from 
Reference  3,  and  the  basic  data  is  tabulated  in  Appendix  II.  The 
vertical  temperature-pressure  profile  at  every  point  along  the  cross 
section  was  assumed  identical  to  that  measured  at  the  MC  location. 

The  lower  elevation  angle  is  limited  to  a  minimum  of  5.5  degrees 
by  the  data  available.  This  can  be  seen  by  sketching  a  ray  at  a  5.5  = 
degree  elevation  angle  starting  at  the  MISTRaM  West  (MW)  site  (Figure  6) ; 
this  ray  reaches  an  altitude  of  10,000  at  MC.  There  is  insufficient  data 
east  of  the  MC  point  for  making  valid  calculations  at  lower  elevation 
angles . 

The  magnitude  and  characteristic  of  the  dispersive  effect  on 
measurements  made  by  the  LIR  are  for  the  environmental  conditions 
measured  at  the  ETR  for  the  3  specific  days.  No  attempt  was  made  to 
determine  the  applicability  of  results  to  other  days  during  the  year. 

The  calculated  values  may  be  high  in  the  case  of  a  few  cloud  formations, 
because  of  low  temperatures  at  high  altitudes,  but  it  is  believed  that 
any  error  resulting  from  the  assumptions,  made  when  valid  data  were  not 
available,  will  provide  an  upper  bound  on  the  magnitude  of  the  effect. 

The  results  tend  to  show  that  the  presence  of  some  clouds  does  not  negate 
the  application  of  this  advanced  technique  for  metric-range  tropospheric 
corrections . 
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Figure  6.  Paths  of  Rays  at  Various  Elevation  Angles,  Cloud  Cross  Section  for  August  7,  1963 


SECTION  III 


COMPUTER  PROGRAMS 

The  computation  problem  was  divided  into  two  phases  because  of 
limitations  on  the  available  amount  of  computer  memory  storage.  The 
first  phase  (LIR)  determines  the  dispersivity  at  each  grid  intersection 
point  in  the  cross  section  and  the  total  dispersivity  for  vertical  rays 
at  points  500  feet  apart  along  the  baseline.  The  flow  diagram  for  LIR 
is  shown  in  Figure  7.  The  dispersivity  distribution  for  vertical  rays 
for  the  3  test  days  has  been  plotted  as  a  function  of  location  along 
the  MISTRAM  baseline  in  Figures  8,  9,  10  and  11.  The  program  is  pre¬ 
sented  in  Appendix  III,  and  tabulated  results  are  presented  in  Tables 
II,  III,  IV,  and  V.  The  total  dispersivity  for  rays  at  elevation  angles 
of  45,  26.6,  18.4,  14.05,  11.3,  9.45,  8.15,  7.1,  6.3,  5.5,  4.8,  4.2, 

4.1,  and  3.8  degrees  is  calculated  in  the  second  phase,  using  the  point 
dispersivity  values  calculated  in  the  LIR  program.  This  second  phase 
is  designated  PHS,  signifying  phase  shifts,  to  avoid  confusion  with 
the  LIR  phase.  The  PHS  program  flow  diagram  is  shown  in  Figure  12. 

In  addition,  provisions  were  made  to  determine  the  total  dispersivity 
along  a  ray  at  a  given  angle  for  any  starting  point  along  the  baseline. 
Movement  of  the  starting  point  of  the  ray  along  the  baseline  is  re¬ 
ferred  to  as  slip  and  is  in  steps  of  500  feet.  The  dispersivity  for 
constant  elevation  angles  of  18.6,  14.0,  and  8.1  degrees  is  plotted  in 
Figures  13  and  14  for  August  7  and  8 , respectively .  The  dashed  portion 
of  the  curve  is  included  to  give  some  additional  indication  of  the 
behavior  of  the  curves,  but  the  basic  data  is  not  adequate  to  insure 
that  this  is  the  total  dispersivity;  the  true  value  may  be  higher  than 
the  value  shown.  The  PHS  program  listing  is  included  in  Appendix  I,  and 
tabulated  values  of  the  dispersivity  as  a  function  of  elevation  angle  for 
10  locations  along  the  baseline  are  included  in  Appendix  II. 
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TABLE  II 


Dispersivity  for  Vertical  Rays  for  Positions  Along  the  MISTRaM  E-W  Base¬ 
line,  August  7,  1963. 


Location 

Degrees 

Location 

Degrees 

Location 

Degrees 

18.0 

0.05 

62.0 

0.17 

94.5 

1.67 

18.5 

0.23 

62.5 

0.40 

95.0 

1.30 

19.0 

0.36 

63.0 

0.48 

95.5 

1.25 

19.5 

0.44 

63.5 

0.48 

96.0 

1.57 

20.0 

0.44 

64.0 

0.48 

96.5 

1.44 

20.5 

0. 36 

64.5 

0.46 

97.0 

0.75 

21.0 

0.47 

65.0 

0.40 

97.5 

0.56 

21.5 

0.48 

65.5 

0.35 

98.0 

0.58 

22.0 

0.49 

66.0 

0.26 

98.5 

0.56 

22.5 

0.32 

66.5 

0.26 

99  .0 

0.78 

23.0 

0.05 

67.0 

0.23 

99.5 

0.78 

27.0 

0.01 

67.5 

0.06 

100.0 

0.87 

27.5 

0.09 

76.0 

0.01 

100.5 

0.70 

28.0 

0.25 

76.5 

0.01 

101.0 

0.45 

28.5 

0.29 

77.0 

0.26 

29.0 

0.30 

77.5 

0.35 

29 .5 

0.30 

78.0 

0.44 

30.0 

0.13 

78.5 

0.47 

30.5 

0.01 

79.0 

0.49 

40.5 

0.01 

79.5 

0.51 

41.0 

0.12 

80.0 

0.52 

41.5 

0.17 

80.5 

0.54 

42.0 

0.16 

81.0 

0.50 

42.5 

0.15 

81.5 

0.25 

43.0 

0.14 

82.0 

0.24 

43.5 

0.18 

82.5 

0.46 

44.0 

0.13 

83.0 

0.51 

45.5 

0.02 

83.5 

0.50 

46.0 

0.06 

84.0 

0.50 

46.5 

0.03 

84.5 

0.49 

55.0 

0.01 

85.0 

0.30 

55.5 

0.06 

85.5 

0.16 

56.0 

0.10 

86.0 

0.01 

56.5 

0.13 

91.0 

0.22 

57.0 

0.12 

91.5 

0.97 

57.5 

0.00 

92.0 

1.55 

58.5 

0.01 

92.5 

1.66 

59.0 

0.12 

93.0 

1.78 

59.5 

0.13 

93.5 

2.03 

K  60.0 

0.06 

94.0 

2.01 

- 
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TABLE  III 


Dispersivity  for  Vertical  Rays  for  Positions  Along  the  MISTRaM  E-W  Base¬ 
line,  August  8,  1963. 


Location 

Degrees 

Location 

Degrees 

Location 

Degrees 

Location 

Degrees 

18.0 

0.52 

45.5 

0.28 

84.0 

1 . 28 

105.0 

4.22 

18.5 

1.00 

46.0 

0 . 66 

84.5 

1.36 

105.5 

5.02 

19.0 

1.18 

46.5 

1.38 

85.0 

1.46 

106.0 

4.48 

19.5 

1.30 

47.0 

2.58 

85.5 

1.58 

106.5 

5.18 

20.0 

1.26 

47.5 

3.78 

86.0 

1.38 

107.0 

5.69 

20.5 

1.32 

48.0 

4.71 

86.5 

0.57 

107.5 

6.03 

21.0 

1.21 

48.5 

5.79 

87.0 

0.24 

108.0 

6.36 

21.5 

0.95 

49.0 

6.63 

87.5 

0.23 

108.5 

6.18 

22.0 

0.53 

49.5 

6.89 

88.0 

0.43 

109.0 

5.48 

25.0 

0.18 

50.0 

7.20 

88.5 

0.66 

109.5 

3.47 

25.0 

0.41 

50.5 

7.23 

89.0 

0.49 

110.0 

0.83 

26.0 

0.52 

51.0 

7.42 

89.5 

0.57 

110.5 

0.65 

26.5 

0.57 

51.5 

7.32 

90.0 

0.90 

111.0 

0.18 

27  .0 

0.54 

52.0 

6.74 

90.5 

2.15 

111.5 

0.01 

27.5 

0.24 

52.5 

6.48 

91.0 

2.29 

112.0 

0.39 

33.5 

0.05 

53.0 

6.30 

91.5 

2.37 

112.5 

0.62 

34.0 

0.18 

53.5 

6.45 

92.0 

2.79 

113.0 

0.69 

34.5 

0.32 

54.0 

5.79 

92.5 

3.19 

113.5 

0.73 

35.0 

0.53 

54.5 

4.42 

93.0 

3.21 

114.0 

0.74 

35.5 

0.88 

55.0 

1.52 

93.5 

3.03 

114.5 

0.69 

36.0 

1.31 

55.5 

0.05 

94.0 

3.27 

115.0 

0.41 

36.5 

1.73 

63.5 

0.02 

94.5 

3. 24 

115.5 

0.06 

37.0 

2.27 

64.0 

0.15 

95.0 

2.76 

116.5 

0.01 

37.5 

2.69 

64.5 

0.34 

95.5 

1.94 

117.0 

0.04 

38.0 

3.03 

65.0 

0.53 

88.5 

0.03 

117.5 

0.37 

38.5 

3.00 

65.5 

0.56 

89.0 

0.09 

118.0 

0.77 

39.0 

3.43 

66.0 

0.53 

96.0 

0.62 

118.5 

0.85 

39.5 

3.78 

66.5 

0.24 

96.5 

0.96 

119.0 

0.69 

40.0 

4.02 

67.0 

0.08 

97.0 

1.07 

119.5 

0.40 

40.5 

3.94 

72.5 

0. 29 

97.5 

1.18 

120.0 

0.43 

41.0 

3.77 

73.0 

0.41 

98.0 

1.09 

120.5 

0.37 

41.5 

2.93 

77.0 

0.12 

98.5 

0.72 

121.0 

0.38 

42.0 

2.47 

77.5 

0.34 

99.0 

0.42 

121 .5 

0.21 

42.5 

3.53 

78.0 

0.10 

99.5 

0.08 

122.0 

0.40 

43.0 

3.46 

81.5 

0.00 

102.5 

0.00 

122.5 

0.28 

43.5 

3.02 

82.0 

0.18 

103.0 

0.09 

44.0 

2.04 

82.5 

0.51 

103.5 

0.24 

44.5 

0.96 

83.0 

0.92 

.  104.0 

1.52 

45.0 

0.32 

83.5 

Lk - 

1.12 

104.5 

2 . 66 
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TABLE  IV 


Dispersivity  for  Vertical  Rays  for  Positions  Along  the  MISTRrtM  E-W  Base¬ 
line,  August  9,  1963. 


Location 

Degrees 

Location 

Degrees 

Location 

Degrees 

Location 

Degrees 

52.5 

0.36 

78.5 

1.13 

20.5 

0. 20 

43.0 

0.52 

53.0 

1.06 

79  .0 

1.06 

21.0 

0.16 

43.5 

0  .05 

53.5 

2. 34 

79.5 

0.96 

21.5 

0.02 

52.0 

0.09 

54.0 

3.51 

80.0 

1  .04 

24.5 

0.06 

97.5 

0.28 

54.5 

6.33 

80.5 

1.04 

25.0 

0.40 

98.0 

0.56 

55.0 

8. 29 

81.0 

0.98 

25.5 

1.02 

98.5 

0.54 

55.5 

11.57 

81.5 

0.90 

26.0 

1.88 

99.0 

0.51 

56.0 

15.20 

82.0 

0.79 

26.5 

3.13 

99.5 

0.49 

56.5 

18.65 

82.5 

0.45 

27.0 

4.77 

100.0 

0.18 

57.0 

20 .43 

83.0 

0.22 

27.5 

6.82 

100.5 

0  .05 

57.5 

22.77 

83.5 

0.17 

28.0 

7.44 

10  3.0 

0 . 14 

58.0 

24.58 

84.0 

0 . 16 

28.5 

7.29 

102.5 

0.10 

58.5 

26.66 

84.5 

0.15 

29  .0 

6.06 

103.0 

0.03 

59.0 

30.40 

85.0 

0.07 

29 .5 

5.28 

103.5 

0.14 

59.5 

32.47 

89.5 

0.08 

30.0 

4.77 

104.0 

0.09 

60  .0 

32.54 

90.0 

0.22 

30.5 

4.56 

118.0 

0.01 

60.5 

30.89 

90.5 

0.16 

31.0 

4.38 

118.5 

0 .12 

61.0 

21.62 

91.0 

0.01 

31.5 

3.42 

119 .0 

0.10 

61.5 

14.98 

96.0 

0.01 

32.0 

3.07 

128.0 

0 . 26 

62.0 

13.28 

96.5 

0.09 

32.5 

2.22 

129.0 

2.19 

62.5 

11.63 

97.0 

0.13 

33.0 

1.41 

129.5 

2.97 

63.0 

12.82 

97.0 

0.18 

33.5 

1.26 

130.0 

2.78 

63.5 

10.88 

8.5 

0.11 

34.0 

1.33 

130.5 

3.05 

64.0 

10.15 

9.0 

0.22 

34.5 

1.51 

131.0 

2.25 

64.5 

11.50 

9.5 

0.26 

35.0 

2.17 

131 .5 

0.07 

65.0 

7.22 

10.0 

0.29 

35.5 

2.55 

135.5 

0.03 

65.5 

3.68 

10.5 

0.30 

36.0 

2.56 

136.0 

0.13 

66.0 

1.35 

11.0 

0.33 

36.5 

2.24 

136.5 

0 . 21 

65.5 

4.43 

11.5 

0.56 

37.0 

1.76 

137.0 

0.39 

66.0 

5.60 

12.0 

0.80 

37.5 

1.31 

137.5 

0.49 

66.5 

6.50 

12.5 

0.89 

38.0 

0.83 

138.0 

0.49 

67.0 

6.44 

13.0 

0.94 

38.5 

0.65 

138.5 

0.61 

67.0 

0.11 

13.5 

0.90 

39.0 

0.28 

139 .0 

0.72 

67.5 

2.94 

14.0 

0.90 

39.5 

0.18 

139 .5 

1 . 20 

68.0 

1.86 

14.5 

0.76 

40.0 

0.34 

140  .0 

1.97 

68.5 

0.15 

15.0 

0.18 

40.5 

0.84 

140.5 

2.51 

69.0 

0.00 

18.5 

0.04 

41.0 

1.01 

141 .0 

4.53 

77.0 

0.29 

19  .0 

0.13 

41.5 

1.11 

141.5 

6.11 

77.5 

0.69 

19.5 

0.22 

42.0 

1.13 

142.0 

7.15 

78.0 

1 .00 

20  .0 

0.24 

42.5 

1.06 

142.5 

7.85 
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TABLE  IV  (Cont’d.) 


Location 

Degrees 

143.0 

8. 26 

143.5 

8.59 

144.0 

5.20 

144.5 

4.29 

145.0 

4.20 

145.5 

1.94 

146.0 

0.66 

146.5 

0.32 

147.0 

1.79 

147.5 

1.75 

148.0 

1.72 

148.5 

1.28 

149.0 

0.74 
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Dispersivity  for  Rays  at  Given  Elevation  Angles 
for  10  Positions  Along  the  MISTRAM  E-W  Baseline 
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SECTION  IV 


DISCUSSION  OF  RESULTS 

Variation  of  the  dispersivity  for  vertical  rays  as  a  function  of 
position  along  the  baseline  for  August  7,  8,  and  9,  is  given  by  the 
dispersivity  distribution  curves  in  Figures  8,  9,  10,  and  11.  The 
magnitude  of  the  effect  for  the  vertical  rays  varies  between  zero 
and  a  maximum  of  about  32  degrees  relative  phase  shift,  or  7.55  feet 
error  in  integrated  ref ract ivity .  No  noticeable  common  pattern  was 
discernible  in  the  distributions  covering  the  3  days  that  indicate 
minimums  occur  at  specific  locations.  Cloud  height  has  a  very 
important  influence  on  the  magnitude  of  the  dispersive  effect  since 
the  temperature  decreases  with  increasing  height.  The  decrease  in 
temperature,  coupled  with  the  strong  temperature  sensitivity  indicated 
by  the  dispersivity- temperature  relationship  given  in  Figure  5,  results 
in  relatively  high  values  of  dispersivity.  For  example,  the  temperature 
was  -25.8  degrees  at  27,000  feet  altitude  on  August  9,  according  to  the 
Rawinsonde  data.  However,  it  was  assumed  that  the  water  was  at  air 
temperature  since  no  information  is  available  on  the  actual  temperature 
of  the  measured  liquid  water  at  this  altitude,  nor  on  the  percentage  of 
water  in  the  form  of  ice  crystals  or  snow.  Since  nonresonant  absorption 
of  the  snow  and  ice  crystals  occurs  in  the  kilohertz  and  low  megahertz 
range,  the  particles  are  not  expected  to  be  dispersive  at  15.6  and 
31.2  gHz.  On  the  other  hand,  liquid-coated  particles  would  have  an 
effect  similar  to  the  effect  of  a  corresponding  water  droplet. 

In  addition,  little  experimental  data  is  available  on  the  char¬ 
acteristics  of  water  at  the  supercooled  temperature.  The  lowest  tem¬ 
perature  at  which  the  water  characteristics  are  available  was  at  -8°C 
(Reference  3,  pages  13  to  17).  Hence,  great  caution  should  be  exercised 
in  using  the  maximum  values  that  occurred  in  the  high  clouds.  The 
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calculated  values  can  be  considered  as  upper  bounds;  if  one-half  the 
measured  liquid  was  in  the  form  of  dry  H^O  ice  crystals,  then  the  mag¬ 
nitude  of  the  effect  should  be  decreased  by  50  percent. 

Variation  of  the  dispersivity  as  a  function  of  elevation  angle  is 
shown  in  the  graphs,  Figures  15  and  16 ,  for  August  7  and  8,  respectively. 
It  is  evident  that  the  dispersivity  is  erratic  and  varies  from  day  to  day 
and  with  location  along  the  baseline. 

The  characteristic  common  to  plots  of  this  data  type  is  elevation 
angles  at  which  the  dispersive  effect  has  very  low  values,  the  mag¬ 
nitude  generally  increasing  with  decreasing  elevation  angles.  For 
example,  the  error  in  the  integrated  refractivity  graph  at  both  the 
11.3-  and  14. 0-degree  elevation  angle  points  on  August  8,  1963  (Figure 
16)  was  about  0.9  foot,  after  reaching  a  maximum  of  2-4  feet  at  25.6 
degrees.  While  no  calculations  were  performed  in  this  investigation  to 
determine  the  magnitude  of  attenuation  of  water  droplets  in  clouds,  it 
is  obvious  from  physical  considerations  that  amplitude  of  the  signals 
would  decrease  in  regions  where  magnitude  of  the  differential  phase 
increases.  Hence,  the  possibility  exists  of  monitoring  amplitude  of  the 
signals  to  obtain  some  indication  of  when  the  minimum  error  in  integrated 
refractivity  occurs.  The  fact  that  the  dispersivity  tends  to  be  larger 
at  the  lower  elevation  angles  is  expected  because  the  probability  of 
encountering  clouds  increases  with  decreasing  elevation  angle.  The 
graphs  of  dispersion- vs .- ray  starting  point,  in  Figures  13  and  14,  on 
the  average  also  tend  to  show  the  same  effect. 

Another  limitation  in  this  study  was  the  lack  of  data  covering  a 
sufficiently  larger  geographical  area  for  various  times  during  one  day. 
The  photographs  of  cloud  structure  (in  Figure  XIV  and  on  page  40, 
Reference  6)  indicate  that  there  were  few  clouds  east  of  MISTRAM 
CENTRAL.  The  calculated  value  of  effective  cloud  signals  for  rays  at 
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Figure  16.  Calculated  Effect  of  Liquid  Water  Droplets  Upon  LIR,  August  8,  1963,  for  Ray  Starting 
at  M1STRAM  West  Location 


elevation  angles  of  less  than  5.5  degrees  would  be  lower  than  the  actual 
value  because  the  values  of  measured  liquid  water  would  appear  as  zero 
for  all  positions  east  of  the  last  measurement  points.  This  same  effect 
occurs  when  the  starting  point  of  the  rays  is  moved  along  the  baseline 
from  the  MW  to  MC  location,  or,  in  other  words,  when  the  slip  value 
increases  to  large  values.  This  is  the  reason  that  the  righthand  ends 
of  the  curves  in  Figures  13  and  14  are  dashed  lines.  For  example,  it  is 
evident  in  the  cloud  cross  section  for  August  7  (Figure  6)  that  the  ray 
at  5.5  degrees,  that  is  N  =  10,  reaches  10,000  feet  at  the  MC  location. 
The  calculated  value  of  the  effect  on  any  ray  at  a  lower  elevation  angle 
would  not  have  the  correct  value  if  there  were  clouds  east  of  MC  and 
below  10,000  feet  altitude  which  were  not  measured. 

The  vast  amount  of  data  for  August  9  resulted  in  computer  storage 
problems  during  the  PHS  program  and  the  results  were  not  included  since 
they  could  be  misleading.  The  calculated  data  is  in  terms  of  degrees  of 
differential  phase  shift  referred  to  31.2  gHz.  Each  degree  of  phase 
shift  referred  to  the  31.2  gHz  frequency  corresponds  to  a  0.236-foot 
error  in  the  distance  measured,  an  error  in  integrated  ref ractivity . 

Numerical  values  of  the  dispersivity ,  as  a  function  of  elevation 
angle  and  position  500  feet  apart  of  the  ray  starting  point  along  the 
baseline,  are  included  in  Tables  II  and  III  for  August  7  and  8.  The 
values  for  elevation  angles  45,  18.6,  14  and  8.1  degrees  are  plotted  in 
Figures  13  and  14  and  show  that  no  consistent  pattern  in  the  relative 
positions  of  the  curves  is  evident. 


SECTION  V 


CONCLUSIONS 

1.  The  presence  of  cloud  droplets  in  the  measurement  path  would  intro¬ 
duce  errors  in  LIR  measurements. 

2.  The  magnitude  of  error  is  a  function  of  the  elevation  angle  and  type 
of  cloud  formation. 

3.  There  are  many  more  so-called  holes  in  the  cloud  structure  at  the 
higher  elevation  angles  than  at  the  lower  elevation  angles;  at  the 
higher  angles  the  error  introduced  by  the  could  droplets  is  small 
compared  to  the  total  refraction  correction.  For  example,  there  is 
a  relative  phase  shift  of  4  degrees,  or  a  refraction  error  of  only 
0.9  foot  (Figure  16  at  14  degrees)  out  of  a  total  refraction  error 
of  approximately  25  feet  for  a  target  at  about  14  degrees  elevation 
angle,  or  at  20-km  height  and  700-km  range  (NBS  7254).  The  residual 
error  is  about  2  feet  RMS  for  this  condition. 

4.  Field-test  experience  indicates  that  signal  amplitude  can  be  used 

as  an  indicator  of  the  presence  of  the  cloud-droplet  error  minimums. 
Hence,  it  may  be  possible  to  monitor  the  data  and  make  corrections 
based  on  the  special  data  points. 

5.  Further  studies  are  necessary  if  it  is  desired  to  determine  the 
applicability  of  this  study  to  other  sites. 


34 


APPENDIX  I 


AN  ESTIMATE  OF  THE  MAGNITUDE 
OF  MICROWAVE  DISPERSION  IN  CLOUDS  AND  RAIN* 


A  normal  atmosphere  with  water  vapor  as  an  ingredient  will  exhibit 
dispersion  in  the  vicinity  of  the  rotational  absorption  line  of  22,235 
me  per  sec.^*^  This  absorption  is  caused  by  rotational  transitions 
between  molecular  states  near  1.5  x  10^  me  per  sec.  Liquid  and  solid 
water  structure  prevents  this  rotation,  and  the  dispersion  in  these  forms 
will  be  of  different  magnitude  and  quality. 


The  wavelengths  of  interest  here  are  1.8  cm  (16.7  kmc  per  sec)  and 
0.9  cm  (33.4  kmc  per  sec).  These  wavelengths  are  large  compared  to 
droplets  in  fair  weather  clouds. Droplets  generally  average  0.0010 
-0.0020  cm  diameter.  Rain  drops  can  be  almost  as  large  as  the  short 
wavelength  but  are  generally  less  than  0.5  cm  because  of  the  breakup 
from  aerodynamic  forces.  As  shown  in  Figure  17^J  the  most  probable 
diameter  for  rain  drops  associated  with  a  fall  rate  of  0.5  inch  per  hour 
is  approximately  0.2  cm,  while  the  diameter  for  a  rain  rate  of  4.0  inches 
per  hour  is  0.3  cm. 


In  order  to  estimate  the  index  of  refraction  for  mixtures  of  drop¬ 
lets  and  drops  of  water  in  air,  it  is  necessary  to  consider  the  polar¬ 
izability  of  the  droplets.  Classical  electrodynamic  theory  as  developed 

[91 

by  Mie  and  others,  Jcirca  1908,  indicates  that  for  spheres,  whose 
diameter  is  small  compared  with  a  wavelength,  the  response  to  electro¬ 
magnetic  fields  is  equivalent  to  an  electric  dipole  of  moment  given 
by 


3  n  -  1 

1 


(8) 


*  This  appendix  was  originally  published  as  The  MITRE  Corporation, 
Bedford,  Mass.,  document  W-5234,  by  H.  M.  Richardson. 
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Figure  17.  Size  Distributions  of  Raindrops  for  Different  Rain  Intensities 


where , 

a  =  the  radius  of  the  sphere, 

n  =  the  complex  index  of  refraction  of  the  material  in  the  droplet, 

and 

2 

n  =  the  complex  dielectric  constant. 


Since  dielectric  properties  of  various  media  are  additive,  it  is  assumed 
that  an  average  dielectric  constant  for  a  cloud  with  N  droplets  per  unit 
volume  can  be  defined  analogously  to  the  definition  for  an  isotropic 
medium,  that  is 


K  -  1 
av 

4tt 


(9) 


as  compared  with 


K.  -  1 

l _ 

4tt 

where,  av  denotes  average,  i 
constant,  E  is  electric  field 
per  unit  volume). 


E  =  P  ,  (10) 

denotes  isotropic,  K  is  dielectric 
— » 

,  and  P  is  polarization  (dipole  moment 


For  N  droplets  per  unit  volume,  the  polarization  is 


P 

a 


-2  ,  ~  _ 

„  n  -  1  3  tT 

N  — -  a  E. 

n  +  2 


(11) 


From  Equations  (9)  and  (11),  and  considering  droplets  where 
is  complex  dielectric  of  the  droplets 


K 
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K  is  near  1,  hence, 
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This  relationship  indicates  that  the  refractivity  of  a  cloud  of  droplets 
depends  essentially  on  the  percent  volume  content  of  water  with  a  cor¬ 
rection  for  the  dielectric  constant  for  liquid  water. 


Equation  (13)  may  be  rewritten  as 


n _ 

av 


M 

P 


(14) 


where , 

M  =  water  density  of  the  dispersion,  and 
p  =  liquid  water  density. 

3  3 

For  M  in  micrograms  per  cm  and  p  in  grams  per  cm  ,  the  second  term 

3 

is  in  refractivity  units,  that  is,  10  micrograms  per  cm  corresponds  to 

_  2 

15  N  (refractivity)  units  for  large  n^ 


The  dielectric  properties  of  water  liquid  are  known  to  vary  near 
1  cm  vacuum  wave  length . L ^ 5  ^ ^  The  variation  of  the  real  part  of  the 
dielectric  function  as  derived  from  References  3  and  10  is  illustrated 
in  Figure  18,  while  Figure  19  shows  the  variation  of  the  imaginary  (or 
loss)  part  of  the  function.  Since  the  wavelength  in  media  with  complex 
dielectric  functions  depends  on  the  real  part  (while  the  attenuation 
by  absorption  in  the  water  depends  on  the  imaginary  part),  the  dis¬ 
persion  at  the  two  frequencies  of  interest  is  estimated  by  the  differ¬ 
ence  between  the  real  parts. 


Dispersion  = 


1 

-) 

2 


1.8  cm 


1 

-) 

2 


a9 


(15) 


where , 


R  = 


Real  part  of 

n,  =  m  -  iK, 
d 
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Figure  18.  Real  Components  of  Index  of  Refraction 
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Figure  19.  Imaginary  Component  of  Index  of  Refraction 


?  2  2 

Kz  +  2)  +  4m  K 


(16) 


R  = 


2  2  2 
(m  -  KZ  -  1)  (m 


2  2  2  2  2 
Cm  -  KZ  +  2)Z  +  4ni  KZ 


Using  values  obtained  by  interpolation  and  extrapolation  of  the 

3 

data  from  the  two  sources,  M  =  10  micrograms  per  cm 

[10] 


Dispersion  =  2.5  at  1.5  degrees 

=0.0  N  at  0  degrees  C 


[3] 


These  values  may  be  compared  with  0.03  N  units  for  vapor  dispersion 

3 

as  determined  in  Reference  7,  for  10  grams  per  meter  and  20  degrees  C. 

The  data  are  considered  more  reliable  since  the  maxima  of  the 

imaginary  component  occurs  nearer  in  frequency  to  the  half -value  point 
of  dielectric  constant  as  indicated  by  relaxation  theory.  Further,  the 
lack  of  dispersion  for  the  data  from  Reference  3  is  not  physically 
reasonable.  The  variation  of  dispersion  with  temperature  for  Ref¬ 
erence  10  data  is  illustrated  in  Figure  20. 

SNOW  AND  ICE  CLOUDS 

Snow  and  ice  particles  are  not  expected  to  be  dispersive  at  the 
frequencies  of  interest  as  the  nonresonant  absorptions  occur  in  the  kc 
and  me  range.  However,  partially  melted  or  liquid-coated  particles 
should  approach  values  for  water  drops. 


CONCLUSION 

This  preliminary  estimate  obviously  does  not  include  all  the  effects 
to  be  expected  by  scattering  the  multiple  scattering  of  electromagnetic 
energy  by  real  distribution  of  drops  in  a  finite  beam  of  microwave 
energy  traversing  a  turbulent  cloud. 

However,  it  is  apparent  that  dispersion  will  exist  to  at  least  the 
degree  estimated  and  that  more  complete  theory  would  indicate  more 
dispersion  of  the  same  sign. 
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I A  - 19,632| 


TEMPERATURE 


Figure  20.  Percent  Dispersion  vs.  Temperature 
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APPENDIX  II 


basic  data 


Liquid  Water  Droplet  Field  Test  Data  -  August  7,  8,  and  9,  1963 
page  54  . 

Rawinsonde  Data  August  7,  8,  and  9,  1963  page  44  . 
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LIQ  WATER  DROPLET  FIELD  TEST  DATA  AUGUST  7.1963 
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